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Abstract: Interference of standing waves in electromagnetic resonators forms the 
basis of many technologies, from telecommunications [1] and spectroscopy [2] to de¬ 
tection of gravitational waves [3]. However, unlike the confinement of light waves 
in vacuum, the interference of electronic waves in solids is complicated by boundary 
properties of the crystal, notably leading to electron guiding by atomic-scale poten¬ 
tials at the edges [4-7]. Understanding the microscopic role of boundaries on co¬ 
herent wave interference is an unresolved question due to the challenge of detecting 
charge flow with submicron resolution. Here we employ Fraunhofer interferometry 
to achieve real-space imaging of cavity modes in a graphene Fabry-Perot (FP) res¬ 
onator, embedded between two superconductors to form a Josephson junction [8]. 
By directly visualizing current flow using Fourier methods [9], our measurements re¬ 
veal surprising redistribution of current on and off resonance. These findings provide 
direct evidence of separate interference conditions for edge and bulk currents and 
reveal the ballistic nature of guided edge states. Beyond equilibrium, our measure¬ 
ments show strong modulation of the multiple Andreev reflection amplitude on an off 
resonance, a direct measure of the gate-tunable change of cavity transparency. These 
results demonstrate that, contrary to the common belief, electron interactions with 
realistic disordered edges facilitate electron wave interference and ballistic transport. 


Graphene provides an appealing platform to explore “electron-optics” due to the ballistic nature of 
wavelike carriers and ability to engineer transmission of electronic waves in real space using electro¬ 
static potentials [10-17]. In particular, the electronic analog to refractive index is the Fermi energy, 
which is tunable via electrostatic gating [11,18]. Because the gapless spectrum of Dirac materials 
enables continuous tunability of carrier polarity, positive and negative index of refraction regions can 
be combined in bipolar structures that form the building blocks of Veselago “electronic lenses” [15], 
Fabry-Perot (FP) interferometers [11-15,17], and whispering gallery mode cavities [19]. Electronic 
analogs to optical interferometers attract attention because relativistic effects such as hyperlensing 
and phase-coherent Klein transmission provide capabilities beyond conventional optics [10-17,20]. 
Here we investigate the simplest analog to an optical interferometer, the electron FP resonator, which 
consists of standing electron waves confined between two reflective interfaces [21,22]. Despite ex¬ 
tensive exploration in the momentum domain, in which Fermi momentum is simply tuned with a gate, 
little information is available about the real-space distribution of current flow due to the challenge of 
imaging current paths with submicron resolution. Furthermore, in real devices, atomically sharp po¬ 
tentials at the edges of graphene can confine electron waves into guided edge modes, in analogy to the 
guiding of light in optical fibers [4-7], as we have demonstrated experimentally in prior work [23]. To 
investigate the nature of these boundary currents, we measure the interference of standing waves in a 
graphene Josephson junction and image the real space distribution of supercurrent flow using Fraun¬ 
hofer interferometry [9]. By visualizing the spatial structure of current-carrying states in the cavity 
using Fourier methods, our measurements disentangle edge from bulk current flow and highlight the 
surprising role of the crystal boundaries. 

In a coherent electron cavity, quantum interference of electron waves replaces classical diffusion 
as a key feature of electronic transport [21,22]. In our system, a pair of superconducting electrodes 
is coupled to a graphene membrane, defining a ballistic cavity between the two graphene-electrode 
interfaces. As the Fermi wavelength in the cavity is tuned with a gate, the quantized energy levels of 
the cavity are moved on and off resonance with the Fermi energy of the superconducting leads, thus 
inducing an oscillatory critical current whose period satisfies the FP interference conditions. Due 
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to the chiral nature of fermions in monolayer and bilayer graphene, at zero magnetic field carrier 
trajectories with an incidence angle 6 and refraction angle angle 9' produce a contribution to FP 
fringes in the single-particle transmission probability of the form 
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where ti ^2 and ri ^2 are the angle-dependent transmission and reflection amplitudes for the two p-n 
junctions. The resulting fringes are dominated by the angles for which both the transmission and 
reflection are reasonably high (the first harmonic of FP fringes is at its brightest when the product of 
transmission and reflection coefficients \t{6)\‘^\r{6)\‘^ takes a maximum value). In general, a spread 
of angles for different trajectories in the bulk gives rise to a spread of the FP oscillation periods, 
somewhat reducing the fringe visibility in the net current. In contrast, no suppression is expected for 
interference fringes due to edge modes, as discussed in detail later. 


We employ proximity induced superconductivity to shed light on the microscopic nature of electron 
interference in a graphene Josephson resonator [24-27]. On a practical level, graphene provides an 
accessible interface for superconducting electrodes because it is purely a surface material, unlike 2D 
electron gases embedded in semiconductor heterostructures. Although graphene is not intrinsically 
superconducting, proximity-induced superconductivity can be mediated by phase coherent Andreev 
reflection at the graphene/superconductor interface. This process features an electron-hole conversion 
by the superconducting pair potential that switches both spin and valley to preserve singlet pairing and 
zero total momentum of the Cooper pair [28]. In this study, we employ gated mesoscopic Josephson 
junctions consisting of bilayer graphene suspended between two superconducting Ti/Al electrodes, as 
well as a graphene device on hBN. The superconductors serve three roles: (1) they create electrostatic 
potentials that confine electron waves, serving as electronic analogs to mirrors (2) superconducting 
interferometry can extract spatial information on how current flows through the system, and (3) be¬ 
yond equilibrium, scattering events between the superconductors and graphene (multiple Andreev 
reflections) depend critically on resonance conditions and reveal how the resonator couples to the 
outside world. 


A schematic of a suspended graphene Josephson junction is provided in Figure la. To access the 
ballistic regime, we developed a new method to isolate the flake from charge disorder in the underly¬ 
ing dielectric by suspending it over the back gate electrode, described in detail in the Supplementary 
Methods. This approach combines the high purity of suspended devices with superconductivity en¬ 
ables creation of ballistic waveguides where the mean free path of electrons exceeds channel length 
L. We note that similar results are also obtained on a gate-defined resonator in monolayer graphene 
encapsulated in hBN, discussed later, which enables a higher degree of electronic control over the 
cavity while preserving sample quality. 

The superconducting leads serve not only as electronic probes but also induce a resonant electron 
cavity in the scaling limit le > L (Fig. lb) [21,22]. The graphene in the immediate vicinity of 
the Ti/Al contact is n-doped by charge transfer [29], forming an intrinsic n-n or n-p junction near 
the interface when the graphene has electron or hole carriers, respectively. We exploit contact in¬ 
duced doping to define the resonator because it is scalable to ultrashort channel lengths, provides 
electrostatic barriers that are sharp compared to the electron wavelength, and is less complex than 
gate-defined methods [30-33]. Analogous to an optical Fabry-Perot cavity, the n-p junctions serve 
as the electronic counterparts to mirrors while the ballistic graphene channel serves as an electron 
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waveguide. The Fermi wavelength Xp of eleetrons in the eavity is direetly tunable with a gate elee- 
trode whieh eontrols the earrier density n. 

Fabry-Perot (FP) resonanees in ballistie junetions arise due to refleetion from p-n junetions formed 
near supereondueting leads when earrier polarity in the graphene region is opposite to the polarity 
of eontaet doping. Figure lb shows a plot of the normal resistanee i?„, obtained by sweeping the 
gate voltage I 4 at a fixed bias exeeeding Ic- We observe well-resolved resistanee oseillations at small 
positive earrier densities (14 < 0) when n-p-n junetion formation is favored and monotonie behavior 
when doping is unipolar. The dips in ii„ eoineide with earrier densities satisfying the eonstruetive 
interferenee eondition 2d = mXp for eleetron waves in a resonator, where d is the effeetive eavity 
length and m is an integer. Sweeping the gate voltage ehanges the Fermi energy in the graphene and 
henee the Fermi wavenumber, given by kp = 2p(Xp = y44r for a 2D Fermi disk with fourfold 
degeneraey. The eorrespondenee to FP interferenee eonditions ean be seen more elearly in Supple¬ 
mentary Figure SI, whieh shows that is periodie as a funetion of 2d/Xp. Reprodueibility of the 
oseillation period is demonstrated in three deviees of 500 nm length (Fig. SI), while shorter june¬ 
tions exhibit larger periods as expeeted. Quantum eonfinement between the eavity “mirrors” gives 
rise to diserete energy levels with spaeing hvp/2d, where vp = hkp/m* is the Fermi veloeity and 
and m* is the effeetive eleetron mass in bilayer graphene. We evaluate this energy seale to be of the 
order 1 meV using the height of FP diamonds, as measured using voltage bias speetroseopy (Fig. SI). 

The interplay between eavity resonanees and supereurrent is evident from a resistanee eolormap 
as a funetion of Ipc and I 4 (Fig. le-d) showing eritieal eurrent oseillations whose period satisfies 
FP inferferenee eondifions, eonsisfenf wifh supereurrenf propagafion via ballisfie eharge earriers [34]. 
As Ai? in fhe eavify is funed wifh fhe gafe, fhe quanfum levels of fhe eavify are moved on or ojf reso- 
nanee wifh fhe Fermi energy of fhe supereonduefing leads, fhus indueing a oseillafing erifieal eurrenf 
periodie in y/n for bilayer graphene. This phenomena is observed in fwo independenf systems: (1) 
suspended bilayer graphene resonators defined by eonfaef-indueed doping (Fig. Id) and (2) a gafe- 
defined resonafor in monolayer graphene on hBN (Fig. le and Supplemenfary Fig. S2), bofh of whieh 
exhibif similar behavior. In fofal, five suspended bilayer deviees are sfudied wifh a lifhographie dis- 
fanee L befween supereonduefing eonfaefs of 350 to 500 nm and eonfaef widfh FF of 1.5 fo 3.2 ^m, in 
addifion fo one gale-defined monolayer deviee wifh eavify dimensions of L = 100 nm and FF = 2.7 
(see Supplemenfary Melhods). Figure le displays erifieal eurrenf modulalions in a gale-defined 
monolayer resonafor whose oseillalions are periodie in n, in agreemenl wifh a monolayer FP model 
for eavify lenglh ~ 100 nm. 

Nexl we employ supereonduefing inlerferomelry as a fool fo spalially resolve oplies-like phenom¬ 
ena assoeialed wifh eleelron waves eonfined wilhin a ballisfie graphene Josephson junelion. Unlike 
experimenls in ID systems [24,25,35], one ean Ihread flux Ihrough fhe junelion and explore fhe rieh 
inlerplay befween magnelie inferferenee effeels and eavify Iransmission. Upon appliealion of a mag- 
nelie field B, a flux penelrales fhe junelion area and induees a supereonduefing phase differenee 
A(j){x) = 2 p^x/^qW parallel to fhe graphene/eonlael inlerfaee, where <l>o = h/2e is fhe flux quan- 
lum, h is Planek’s eonslanl, and e is fhe elemenlary eharge. When a flux penelrales fhe junelion area, 
fhe erifieal eurrenf Ic{B) exhibils oseillalions in magnelie field given by: 


Ic{B) 



J{x) ■ 


( 2 ) 
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where L is the distanee between supereondueting eleetrodes (Fig. 1) [9,36]. This integral expression 
applies in the wide junetion limit, relevant for our system, where L W and the eurrent density 
is only a funetion of one eoordinate. Beeause the eritieal eurrent Ic{B) equals magnitude of the 
eomplex Fourier transform of the real-spaee supereurrent distribution J{x), the shape of the interfer- 
enee pattern is determined direetly by the spatial distribution of supereurrent aeross the sample [9,37]. 

To visualize eurrent flow assoeiated with interfering eleetron waves in graphene, we measure su¬ 
pereurrent modulations in B field that arise from a Fraunhofer diffraetion. Figure 2a is a eolor map of 
eritieal eurrent Ic as a funetion of gate voltage and magnetie field. Eaeh pixel is obtained by measur¬ 
ing the DC voltage Vsd aeross the junetion as a funetion as a funetion of applied DC eurrent bias Inc 
and extraeting the maximum derivative dVgd/dloc- In n eonventional graphene Josephson junetion 
with uniform eurrent density, the normalized eritieal eurrent Ic{B) //c(0) = | sin(7r$ /^o)/ (tt'F/< ho) I 
is deseribed by Fraunhofer diffraetion and should be independent of gate voltage. Our results exhibit 
a striking departure from this pieture and feature nodes in Ic{B) /lc{0) as a funetion of both Vj, and B 
(as shown in Supplementary Fig. S3). Figure 2a and Fig. S3 display the different behavior of Ic ver¬ 
sus B on and off resonanee, where the red and green dotted lines indieate gate voltages eorresponding 
to on and off resonanee eonditions, respeetively. (Reprodueibility of this phenomenon in additional 
samples is shown in Supplementary Fig. S4.) Using Eq.(2), one ean extraet an effeetive spatial distri¬ 
bution of the supereurrent J{x) by taking the inverse Eourier transform of the above Ic{B) line plots 
with the teehnique of Dynes and Eulton [9] (see Supplementary Methods). As revealed in Eig. 2b, the 
normalized spatial distribution features bulk-dominated eurrent flow on resonanee and an enhaneed 
edge eurrent eontribution off resonanee. 

Inspired by the relation between the spatial eurrent distribution J(x) and eritieal eurrent Ic{B) in 
Eq. (2), we direetly model the spatial distribution of eurrent paths for bilayer graphene in the EP 
regime (Eig. 2e-e). These ealeulations take into aeeount guided edge modes due to band-bending at 
the erystal boundaries, whieh have been experimentally observed in Ref. [23]. This eleetron guiding 
effeet ean be quantified by an edge potential, whieh is eapable of eonfining earriers to edge-defined 
‘waveguides’ in analogy to the eonfinement of photons in fiber optie eables. Energies of these edge 
states lie outside the bulk eontinuum (Eig. 2e), whieh ensures an evaneseent-wave deeay of earrier 
states into the bulk. The resulting states are effeetively one-dimensional, propagating as plane waves 
along the graphene edges. Applying the EP quantization eondition in the p-n-p region leads to a 
sequenee of EP maxima positioned at kn = 'KnjL, where n is an integer and L represents distanee 
between supereondueting eontaets. These quasi-lD states guided along the edge feature head-on 
transmission and refleetion and henee should produee mueh stronger EP fringes than the bulk states. 

As shown in the theoretieal dispersion in Pig. 2e, the interferenee eonditions in the bulk and at 
the edge should not eoineide due to the differenee in the earrier dispersion at the edge and in the 
bulk as well as due to the angle-dependenee of the EP period for the latter earriers. Henee a gradual 
inerease of doping will trigger repeated switehing between the bulk-dominated and edge-dominated 
regimes, with the eurrent distribution switehing from an approximately uniform to edge like, ae- 
eordingly. Qualitatively, this would be manifested in the dependenee of measured eritieal eurrent on 
applied magnetie field, switehing between Praunhofer and more SQUID-like behavior (Pig. 2a). 

To quantify these phenomena, we model EP resonanees using the approaeh deseribed in Ref. [23]. 
Assuming that the edge potential is suffieiently short-ranged, we approximate it with a delta funetion. 
We obtain the density of persistent eurrent along the edge (ehosen to be along y axis) from the exaet 
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Green’s funetion G in a mixed eoordinate-momentum representation: 

(3) ^ \G{e,x,x\ky)Jy 

where Jy is the operator for eurrent along the edge (see Eq. (8-9) in supplement of Ref. [23]). The 
sum in Eq.(3) runs over the values kn = im/L for one sign of n. Eaeh term in Eq.(3) has poles 
eorresponding to bounded states for the momentum value kn, eaeh of whieh eorresponds a to eurrent 
maximum at the edge (see Eig. 2e). The predieted spatially resolved eurrent density j{x) aeross the 
sample as funetion of energy is shown in Eigure 2e. To translate this into an experimentally observ¬ 
able quantity, we model Eraunhofer interferenee pattern Ic{B) using the theoretieal amplitude and 
spatial distribution of edge modes (see Supplementary Materials). This result, plotted in Eig. 2d, 
eaptures the key features of the data, namely the redistribution of eurrent on and off resonanee as well 
as the suppression of side lobes’ intensity on resonanee. Thus, the measurements are eonsistent with 
a model that features separate EP interferenee of guided-wave edge eurrents, in parallel to interfer¬ 
enee of bulk modes. This further suggests that the quasi-lD edge eurrents previously observed [23] 
have ballistie eharaeter. Despite its simplified nature, whieh negleets disorder and finite tempera¬ 
ture effeets, our model eaptures the essential features of the measurements. While the edge potential 
featured in this simulation aeeommodates a single edge ehannel, we note that the number of guided 
modes may exeeed one for stronger potentials. In this ease, eaeh mode would contribute indepen¬ 
dently to the interference pattern, giving rise to fringes with complicated multi-period structure at the 
edge. 

We employ yet another property of superconductor-normal-superconductor (SNS) systems to gain 
insight into the coupling between the cavity modes with the superconducting reservoirs. Because 
the phenomenon of multiple Andreev reflection (MAR) is known to be extremely sensitive to the 
coupling between electrons in the normal metal and the superconductor, we use voltage bias spec¬ 
troscopy to map out the interplay between MAR oscillation amplitude and cavity transmission (Eig. 
3a,b). The millielectronvolt energy scale associated with EP interference substantially exceeds the 
A1 superconducting gap A, allowing one to study the system close to equilibrium conditions for the 
resonator. A colormap of resistance as a function of applied voltage bias Vdc and gate voltage I 4 
shows modulations due to EP interference (Pig. 3b). Well defined MAR peaks appear at 2A, A, and 
2A/3 when the density is tuned off resonance, while MAR is completely suppressed on resonance, 
as visible in line cuts of resistance on and off resonance in Pig. 3c (additional data sets are provided 
in Supplementary Pigs. S5-S6). It is notable that the amplitude of the multiple Andreev reflections 
depends strongly on cavity resonance conditions, thereby providing a direct measure of the tunable 
coupling between the resonator and the outside world. 

The change in visibility of MAR on and off EP resonances is most naturally explained by changes 
in the distribution of transmission eigenvalues, which can be understood using the following model. 
Because supercurrent is predominately transmitted by bulk modes, as indicated by the Eraunhofer 
interferometry data (Pig. 2), we simplify our analysis by focusing on resonances of bulk states. The 
magnitude of multiple Andreev reflection peaks is small for modes with high transmission probability 
due to the absent suppression of higher order scattering processes [38]. In a EP cavity a larger fraction 
of the current is carried by highly transmitting modes when the cavity is tuned to the resonant wave 
length. In a short junction different modes contribute independently to the current, producing the 
observed multiple Andreev reflection pattern. In our junction, ^ = hvp/A 450 — 700 nm, while 
the junction size is Ri 350 nm, so we expect the short junction limit to qualitatively hold. 
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In order to compare this model to the observed experimental data, we have modeled the current 
through the junction as a sum of contributions of the modes with high (~ 0.9), medium (~ 0.6), 
and low transmission (~ 0.3) coefficient (Fig. S7). This separation was chosen to avoid overfitting, 
while keeping the qualitative features of I-V relationships with different transparencies. We then 
approximate 

(4) I{V) = Y,PiTn)-I{V,Tn), 

n 

with Tn the transmission probability in various channels, p the density of transmission eigenvalues, 
and I{V,T) the contribution of a single mode with transmission probability T to the total current, 
calculated in the short junction limit following Ref. [38]. Fitting the model to the measured conduc¬ 
tance curves on and off resonance (Fig. 3c,d and Supplementary Fig. S7) shows that the junction 
transparency is increased on resonance and suggests good qualitative agreement between this theo¬ 
retical interpretation and the experiment. 

We obtain the estimated contributions of each T„ by fitting the measured I-V traces using the 
Eq. (4) constrained by the condition p{Tn) > 0. The fits show no systematic error, and increas¬ 
ing the number of T„ leads to noisier fits, indicating overfitting. The coefficient p corresponding to 
large transmissions increase, while the ones corresponding to low transmissions decrease whenever 
the system is on resonance, at values of the back gate voltage where the normal state conductance is 
also peaked (see Supplementary Fig. S7). The normal state conductance estimated using our model 
Gn = 90 '^nP{Tn) is Smaller than the measured one for all back gate voltages, which may be due 
to deviations from the short junction theory, or the nonlinear behavior of the p-n junctions. 

In summary, we utilize different aspects of proximity-induced superconductivity, particularly Fraun¬ 
hofer interferometry and Andreev scattering, as new tools to resolve optics-like phenomena associated 
with electron waves confined wifhin a ballistic graphene Josephson juncfion. This enables real-space 
visualization of cavity modes in a graphene FP resonafor, which reveals surprising redisfribufion of 
currenf on and off resonance and provides direcf evidence of fhe ballisfic nafure of guided edge cur- 
renfs. These resulfs consfifule a sfrong deparfure from conventional Josephson behavior in graphene 
and mofivafe furlher explorafion of new effecls af fhe infersecfion of superconducfivily and opfics-like 
phenomena. 
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Eigure legends 

Figure 1. Interplay between superconductivity and the Fabry-Perot interference in a bal¬ 
listic graphene Josephson junction, (a) Gated mesoscopic Josephson junction consisting of bilayer 
graphene suspended between two superconducting Ti/Al electrodes. L is the lithographic distance 
between contacts and W is the junction width. In the presence of magnetic field, a flux fhreads fhe 
juncfion area. A currenf bias is applied befween fhe elecfrodes and fhe volfage drop across fhe de¬ 
vice is recorded. A volfage applied fo fhe back gale eleclrode I 4 lunes fhe Fermi wavelenglh Xp 
in fhe cavily. (b) Plol of fhe normal resisfance, oblained by sweeping fhe gale volfage I 4 al a fixed 
bias exceeding Ic- Dala sels in panels (b-d) are from device Bl. Left inset: Charge Iransfer al fhe 
boundaries of fhe superconducling elecfrodes leads lo inlrinsic n-doped regions near fhe conlacls, 
forming an eleclronic resonalor when fhe bulk is luned lo hole doping. Dips in resisfance appear 
when conslruclive inlerferences conditions in fhe cavily are salisfied, 2L = mXp. Right inset: When 
fhe bulk is luned lo eleclron doping, slanding waves are nol formed, leading lo monolonic resisfance. 

(c, d) Plols of resisfance as a funclion of DC currenf bias and back gale volfage. The critical currenf 
Ic oscillates wilh a period lhal salisfies fhe Fabry-Perol (FP) inlerference condilions, consislenl wilh 
supercurrenl propagalion via ballistic charge carriers. (e)Differenlial resisfance of a gale-defined FP 
resonalor in monolayer graphene on hBN (device Ml), as a funclion of lop gate volfage and DC bias 
currenf when fhe back gate volfage is held fixed al -1.75 V. The critical currenf, defined by fhe widlh 
of zero resisfance region along fhe currenf axis, oscillates wilh fhe same periodicily as normal slale 
resisfance, in agreemenl wilh a FP model for cavily lenglh ~ 100 nm. 

Figure 2. Spatially resolved supercurrent imaging in a ballistic graphene cavity, (a) Plot of 
critical current Ic as a function of back gate voltage I 4 and applied magnetic field B. Each pixel 
was obtained by measuring the DC voltage Vgd across the junction as a function as a function of DC 
current bias Ipc and extracting the maximum derivative dYgd/dlpc- Red and green dotted lines 
indicate on and off resonance conditions, respectively. Data was collected from device B2. (b) Real- 
space normalized supercurrent density distribution J{x)/ Jmaxix) extracted from the Ic{B) data in 
(a) using Fourier techniques (see Supplement for details), (c) Spectrum of bilayer graphene with small 
edge potential, for which one edge mode dominates, (d) Theoretical plot of critical current Ic as a 
function of barrier energy and applied magnetic field in presence of edge modes. Bulk and edge cur¬ 
rents produce distinct FP patterns due to different dispersion laws and angle dependent transmission 
of bulk modes, (e) Theoretical calculation of spatially resolved current density across the sample as 
function of energy. Herepo = Xm*/2h, Eq = pQ(2m* and xq = h/po with m*=0.04 nie (BEG band 
mass) and delta function potential strength A = 0.5 eV-nm (see Ref. [23]). Energies corresponding 
to quantized momenta are represented by horizontal red lines. 

Figure 3. Interplay between multiple Andreev reflections and cavity transmission, (a) Schematic 
illustration of the mechanism of multiple Andreev reflection in a graphene Josephson junction for 
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voltage bias eV = 2A/3, (b) Right panel: A colormap of resistance Rn as a function of applied 
voltage bias Vdc and gate voltage 14 shows modulations due to Fabry-Perot interference. Left panel: 
Derivative plot dRu/dVoc for the data on the right. Data sets in panels (c,d) are from device B3. 
(c) Line cuts of resistance versus DC voltage bias on (14=0.3 V, red curve) and off (14=0.14 V, 
blue curve) resonance. Well defined MAR peaks appear at 2A, A, and 2A/3 when the density is 
tuned off resonance, while MAR is completely suppressed on resonance. (d)Theoretically obtained 
conductance profiles in the short junction limit, as a function of applied bias voltage. The curve corre¬ 
sponding to high transmission, Ghigh (red) is computed for a single mode with transmission 0.9. The 
low transmission curve (blue) is obtained for 4 modes with transmission 0.6. Lower transparencies 
lead to the formation of conductance resonances at bias voltages corresponding to 2 A/3, A, and 2A. 
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Figure 3 
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Supplementary Materials and Methods 

Fabrication and design of ballistic graphene Josephson junctions 


Suspended Josephson junctions'. We investigate suspended Josephson junctions of two types. 

The first type, which corresponds to data shown in Fig. Ib-d and Supplementary Fig. SI, features 
superconducting electrodes in the interior of the flake. Graphene is mechanically exfoliated directly 
onto on a 300 nm Si02 dielectric layer that coats a doped silicon wafer serving as a global back 
gate. Next, thin Cr/Au leads are defined using e-beam lithography in a pseudo-four probe geometry 
in order to make electrical contact to the bilayer graphene device. These contacts are spaced roughly 
1-1.5 /xm apart in order to leave room for the superconducting contacts that will eventually define the 
Josephson junction itself. These Cr/Au (3/30 nm) contacts are deposited using thermal evaporation, 
followed by immersion in acetone for metal liftoff. Next, thick gold electrodes are defined in a way 
that overlaps the outer edges of the thin contacts, thus maintaining electrical contact to the flake. The 
thick electrodes serve a dual purpose: (1) to provide structural support and mechanically hold up the 
entire suspended graphene Josephson junction and (2) provide an electrical connection between the 
Josephson junction and the bondpads. After an evaporation mask is defined using e-beam lithogra¬ 
phy, Cr/Au (3/200 nm) is deposited. To define the Josephson junction, a pair of rectangular Ti/Al 
superconducting contacts are patterned in the interior of the flake and extending over the thin Cr/Au 
leads to maintain electrical contact to the bondpads. The superconducting electrodes are patterned 
using e-beam lithography, followed by thermal evaporation of a 10 nm Ti adhesion layer and a 70 nm 
superconducting A1 layer. Finally, in order to protect the superconductor from degradation in acid 
during the suspension process, a PMMA polymer etch mask is defined over the superconducting con¬ 
tacts using e-beam lithography. After development, the entire chip is immersed in a buffered oxide 
wet etchant to remove 150 nm of the underlying Si02 dielectric layer, leaving the Josephson junction 
fully suspended. Immediately following the etch, the substrate is immersed in methanol, followed by 
an acetone soak to dissolve the PMMA mask, after which the chip is again immersed in methanol and 
dried in a critical point dryer. 

The second type of suspended Josephson junction, which corresponds to the data in Fig. 2-3, 
features superconducting electrodes that extend over the full width of the flake. This Josephson 
geometry is preferable for imaging current flow due to the uniform distance between contacts and 
rectangular junction dimensions. Devices are fabricated on a 300 nm Si02 dielectric layer that coats 
a doped silicon wafer that serves as a global back gate. Bilayer graphene flakes are deposited over 
predefined narrow trenches that are etched into the Si02 with a depth of 150 nm. Next, thin Cr/Au 
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contacts and bondpads are defined using e-beam lithography in a pseudo-four probe geometry in or¬ 
der to make electrical contact to the bilayer graphene device. These contacts are spaced far apart 
on either side of the etched trench in order to leave room for the superconducting contacts that will 
eventually define fhe Josephson juncfion ifself. These Cr/Au (3/30 nm) confacfs are deposifed using 
fhermal evaporafion. The devices are fhen immersed in acefone for mefal liffoff, fransferred im- 
mediafely info mefhanol, and carefully dried using a critical poinf dryer due fhe delicafe nafure of 
suspended graphene membranes. To consfrucf fhe Josephson juncfion, superconducting Ti/Al con¬ 
facfs are pafferned along fhe french edges using e-beam lifhography and wifh widfh large enough fo 
achieve elecfrical confacf wifh fhe Cr/Au leads. The superconducfing confacfs are deposifed using 
fhermal evaporafion wifh fhe following procedure: a 10 nm Ti adhesion layer is deposifed, followed 
by a 50 nm layer of A1 superconducfor. As wifh fhe previous sfep, mefal liftoff is conducted by im¬ 
mersion in acefone and mefhanol, followed by drying in a crifical poinf dryer. The mofivafion for 
using Cr/Au bondpads is fo achieve fhe besf possible elecfrical connecfion fo fhe gold bonding wires 
and sample holder pins. Aluminum, by confrasl, oxidizes upon exposure to air and forms infermefal- 
lic compounds af fhe interface wifh gold bonding wire, which would be expected fo degrade elecfrical 
confacf. Devices are currenf annealed in vacuum af dilufion refrigeralor femperafures in order fo re¬ 
move organic processing residues and enhance qualify. All low femperafure dafa is collected using 
sfandard lockin measuremenf fechniques in a Leiden Cryogenics Model Minikelvin 126-TOF dilufion 
refrigerator wifh a base femperafure of ~ 10 mK. 

Suspended Josephson junction device dimensions: 

Sample Bl: Fig. Ib-d, Supplemenfary Fig. SI: main panel plof, righf insefs; blue curve in fhe lefl 
insef. Type 1 geomefry. Disfance befween superconducfing elecfrodes: 500 nm. Widfh of supercon¬ 
ducfing confacfs (defines fransverse dimension of juncfion): 1.7 pm. 

Sample B2: Fig. 2a-b, Supplemenfary Fig. S3: Type 2 geomefry. Disfance befween superconducf¬ 
ing elecfrodes: 350 nm. Juncfion widfh: 1.7 pm. 

Sample B3: Fig. 3, Supplemenfary Fig. S6: Type 2 geomefry. Disfance befween superconducfing 
elecfrodes: 350 nm. Juncfion widfh: 1.7 pm. Note: Dafa sefs B2 and B3 are from fhe same physical 
device buf are collecfed afler differenl currenf annealing iferafions and fhus have differenl disorder 
configurations. 

Sample B4: Supplemenfary Fig. SI: green curve in fhe leff insef. Type 1 geomefry. Disfance 
befween superconducfing elecfrodes: 500 nm. Widfh of superconducfing confacfs (defines fransverse 
dimension of juncfion): 3.2 pm. 

Sample B5: Supplemenfary Fig. SI: red curve in fhe leff insef. Type 1 geomefry. Disfance befween 
superconducfing elecfrodes: 500 nm. Widfh of superconducfing confacfs: 1.65 pm. 

Sample B6: Supplemenfary Fig. S4-S5: Type 2 geomefry. Disfance befween superconducfing 
elecfrodes: 350 nm. Juncfion widfh: 1.5 ^m. 

Sample B7: Supplemenfary Fig. S4: Type 2 geomefry. Disfance befween superconducfing elec¬ 
frodes: 350 nm. Juncfion widfh: 1.5 pm. Note: Dafa sefs B6 and B7 are from fhe same physical 
device buf are collecfed affer differenf currenf annealing iferafions and fhus have differenf disorder 
configurations. 

Josephson junctions on hBN: To invesfigafe a separafe device design, fhe gate-defined FP cavify, 
we also consider one dual-gated monolayer graphene Josephson juncfion encapsulated in hexangonal 
boron nifride (hBN). By isolating fhe graphene from fhe surface roughness and charge disorder asso- 
ciafed wifh fhe underlying silicon dioxide gale dieleclric, hBN subslrales enable high device qualify 
fo be achieved, which is a crucial ingredienf for observing ballistic charge fransporf. This Josephson 
juncfion has a disfance of 750 nm befween fhe superconducfing confacfs and a flake widfh of 2.7 pm. 
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The top gate length, whieh defines the size of the FP resonator, is ~100 nm. The supereondueting 
eleetrodes eonsist of an adhesion layer of Ti (10 nm) and a supereondueting layer of A1 (60 nm). The 
top gate eonsists of Ti/Au (5/50 nm). The thieknesses of the top and bottom hBN flakes that eneap- 
sulate the graphene are ~19 nm and ~30 nm, respeetively, as measured by atomie foree mieroseopy 
(AFM). 

Encapsulated Josephson junction device dimensions: 

Sample Mi: Fig. le. Supplementary Fig. S2: Distanee between supereondueting eleetrodes: 750 
nm. Junetion width: 2.7 pm. The top gate length: ~100 nm. 

Fourier method for extraetion of supereurrent density distribution 


In order to disentangle edge from bulk eurrent flow through the resonator, we employ the Fourier 
teehniques of Dynes and Fulton to reeonstruet the real-spaee supereurrent distribution from the mag- 
netie interferenee pattern Ic{B). This proeedure, described in detail in Ref. (37), is briefly summa¬ 
rized here. When a magnetic field B is applied perpendicular to the junction area, the critical current 
Ic{B) through a Josephson junction is: 


( 1 ) 


Ic{B) = \Ic{B)\ = 


J{x) exp(27ri(L + lAi)Bx/^o)dx 


where x is the dimension along the width of the superconducting contacts (labeled in Fig. 1), L is the 
distance between contacts, l^i is the magnetic penetration length scale (determined by the London 
penetration depth of the superconductor and flux focusing), and = h/2e is the flux quantum. This 
integral expression applies in the narrow junction limit where L <C W, relevant for our system. 


Observing that Ic{B) represents the complex Fourier transform of the current density distribution 
J(x), one can apply Fourier methods to extract the spatial structure of current-carrying electronic 
states. Because the antisymmetric component of J{x) vanishes in the middle of the junction, the 
relevant quantity for analyzing edge versus bulk behavior is the symmetric component of distribution. 
By reversing the sign of Ic{B) for alternating lobes of the superconducting interference patterns, we 
reconstruct Ic{B) from the recorded critical current. One can determine the real-space current density 
distribution across the sample by computing the inverse Fourier transform: 


( 2 ) 


Jsix) 


Tc{B) exp(27ri(L -|- lAi)Bx/^o)dB 


We employ a raised cosine filter to taper the window at the endpoints of the scan in order to reduce 
convolution artifacts due to the finite scan range Bmin < B < Bmax- This the explicit expression 
used is: 

r Bjnax 

(3) Js{x) ^ / Zc{B) cos'^{7rB/2LB) ex.p{2TTi{L + lAi)Bx/^o)dB 

Bmin 

where n = 0.5 — 1 and = {Bmax — Bmin)/"^ is the magnetic field range of the scan. 



4 


Supplementary Figures 



Back gate voltage (V) 


Figure S 1. Characterization of Fabry-Perot interference conditions in ballis¬ 
tic graphene Josephson junctions. Main panel: Plot of the normal resistance Rn, 
obtained by sweeping the gate voltage 14 at a fixed bias exceeding Ic (repeated here 
from Fig. lb for reference to the other panels). Data was collected from sample 
Bl. Right inset: Fabry-Perot diamonds obtained using voltage bias spectroscopy, 
as shown in color maps of R{Q) and its derivative dRn/dVb, as function of back 
gate voltage 14 and voltage bias Vdc- Data from sample Bl. Left inset: Rn plotted 
versus 2d/\F, where d is the effective junction length and Xp is the Fermi wave¬ 
length. By comparing the junction length L to the effective size d extracted from fits, 
we determined that the contact-doped regions extend at most 100 nm into the chan¬ 
nel, consistent with the results of scanning photocurrent studies. Reproducibility of 
the oscillation period is demonstrated in three devices of length L = 500 nm. The 
blue resistance curve is from sample Bl, the green curve is from sample B4, and the 
red curve is from sample B5 and offset by -250 D. Resonances marked by dips in 
resistance appear when constructive interferences conditions are satisfied. 
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Top gate voltage (V) 


Figure S2. Ballistic resistance oscillations in a gate-defined Fabry-Perot inter¬ 
ferometer on hBN. Normal state resistance of a monolayer graphene device on 
hBN as a function of top gate when the back gate is held fixed at -1.75 V (resistance 
line cut corresponds to the white dotted line in the inset). Data was collected from 
sample Ml. The oscillation period agrees with a Fabry-Perot model with a cavity 
length ~ 100 nm. Inset shows that oscillation occurs in p-n-p and n-p-n regions, 
which is characteristic of Klein tunneling in monolayer graphene. The oscillations 
also suggest the ballistic nature of electronic transport in the locally gated region. 





Back gate voltage (V) 
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Figure S 3 . Nontrivial current flow through a graphene Fahry-Perot resonator, 
as revealed hy Fraunhofer interferometry, (a) Unsaturated map of the critical cur¬ 
rent Ic{B) data from Fig. 2a, plotted over a full color scale range, (h) Plot of nor¬ 
malized critical current Ic{B)/Ic{B = 0) from the data in Fig. 2a, indicating a non¬ 
trivial dependence of Fraunhofer interference on cavity resonances. Red and green 
dotted lines indicate on and off resonance conditions for the cavity, respectively, (c) 
Real-space supercurrent density distribution J{x) extracted from the Fraunhofer in¬ 
terference Ic{B) data in Fig. 2a using Fourier techniques. (Fig. 2b in the main text 
is a plot of the real-space normalized supercurrent density distribution for this data 
set.) Data was collected from sample B2. 



Magnetic field (mT) Magnetic field (mT) 


Figure S4. Dependence of normalized Fraunhofer interference on cavity reso¬ 
nances in additional samples. Plot of normalized critical current Ic{B)/Ic{B = 0), 
indicating nontrivial dependence of Fraunhofer interference on cavity resonances. 
Red dotted lines indicate on resonance conditions for the cavity. Data was collected 
from samples B6 (panel (a)) and B7 (panel (h)), which exhibit qualitatively similar 
behavior to sample B2 in Supplementary Fig. S3. 
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Figure S5. Interplay between multiple Andreev reflection and cavity reso¬ 
nances in an additional device. Experimental resistance profiles obtained using 
voltage bias spectroscopy, measured at fixed back gate voltages 14 tuned to on or off 
resonance conditions for the cavity (red or blue curves, respectively). The amplitude 
of the multiple Andreev reflections, manifested in resistance dips at 2 A /n for integer 
re, is strongly modulated by cavity transmission and thus exhibits suppression when 
carrier density is tuned on-resonance, (a) Red curve: 14 = —0.7 V, (on-resonance, 
corresponding to a dip in normal state resistance). Blue curve: 14 = —0.85 V, 
(off-resonance, corresponding to a peak in normal state resistance), (b) Red curve: 
14 = —0.95 V (on-resonance); Blue curve: 14 = —1.025 V (off-resonance). In 
panels (a) and (b), the red curves are offset by —50011 for clarity, (c) All curves 
from panels (a) and (b), plotted on the same resistance scale. Red curve: 14 = —0.7 
V (on-resonance); dark blue curve: 14 = —0.85 V (off-resonance); magenta curve: 
14 = —0.95 V (on-resonance).; light blue curve: 14 = —1.025 V (off-resonance). 
Data was collected from sample B6, which shows qualitatively equivalent behavior 
to sample B3 in Fig. 3 of the main text. 
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Figure S6. Additional multiple Andreev reflection data sets from sample B3. 

More voltage bias spectroscopy data from sample in Fig. 3, taken at additional back 
gate voltages. Red curve: 14 = —0.1 V (on-resonance); Blue curve: 14 = —0.3 V 
(off-resonance). 
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Figure S7. Theoretical dependence of multiple Andreev reflection amplitude 
on cavity transmission, (a) Theoretical resistance profiles as a function of applied 
bias voltage, corresponding to the experimental data in Fig. 3c. The red curve (on 
resonance, 14 = 0.317) shows suppressed MAR features, while in the blue curve (off 
resonance, I 4 = 0.517) MAR peaks appear at bias voltages 2A, A, and 2A/3. The 
red curve has been shifted upwards by 0.015i?o for clarity, (b) Simulated resistance 
map obtained by fitting the measured data to the short junction model, plotted con¬ 
tinuously as a function of applied DC bias voltage 17 and back gate voltage. The 
theoretical resistance profile is in good agreemenf wifh fhe experimenfal one (Fig. 
3b), showing well defined MAR peaks when fhe sysfem is off resonance, and sup¬ 
pressed MAR feafures on resonance, (c) Mode confribufions p (fhick solid lines) 
corresponding fo large (red), medium (green), and small (blue) fransmissions, as a 
funcfion of back gale vollage. The black and blue dashed lines show fhe normal sfafe 
conducfance values from measuremenf and fheory, respecfively. 













